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The dynamics of particle, event, and apparent horizons in FLRW space are discussed. The 
apparent horizon is trapping when the Ricci curvature is positive. This simple criterion coincides 
with the condition for the Kodama-Hayward apparent horizon temperature to be positive, and also 
discriminates between timelike and spacelike character of the apparent horizon. We discuss also the 
entropy of apparent cosmological horizons in extended theories of gravity and we use the generalized 
2nd law to discard an exact solution of Brans-Dicke gravity as unphysical. 
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INTRODUCTION 

Black hole thermodynamics [l[ links classical gravity 
and quantum mechanics and constitutes a major ad- 
vancement of the theoretical physics of the 1970's. The 
discovery by Bekenstein Q and Hawking d, 0] that black 
hole horizons have entropy and temperature associated 
with them allowed for the formulation of a complete ther- 
modynamics of black holes. It is widely believed that 
formulating also a statistical mechanics to explain black 
hole thermodynamics in terms of microscopic degrees of 
freedom requires a fully developed theory of quantum 
gravity, which is not yet available. 

Soon after the discovery of Hawking radiation 0, Q] , 
Gibbons and Hawking discovered that also the de Sitter 
cosmological event horizon is endowed with a tempera- 
ture and an entropy, similar to the Schwarzschild hori- 
zon Later, it was realized that the notion of black 
hole event horizon, which requires one to know the en- 
tire future development and causal structure of space- 
time, is essentially useless for practical purposes. The 
teleological event horizon is not an easy quantity to com- 
pute: this feature has been emphasized by the develop- 
ment of numerical relativity. In the sophisticated simu- 
lations of black hole collapse available nowadays, outer- 
most marginally trappes surfaces and apparent horizons 
are used as proxies for event horizons [6] . While the early 
literature on black holes and the development of black 
hole thermodynamics in the 1970's focused on static and 
stationary black holes, for which apparent and event hori- 
zons coincide, dynamical situations such as the interme- 
diate stages of black hole collapse, black hole evaporation 
backreacting on its source, and black holes interacting 
with non-trivial environments (e.g., with another black 
hole or compact object, or with a cosmological back- 
ground) require the generalization of the concept of event 
horizon to situations in which no timelike Killing vector 
is available. For this purpose, the concepts of apparent, 
trapping, isolated, dynamical, and slowly evolving hori- 
zons were developed (see [7H9j for reviews). 

In addition to black hole horizons, also cosmological 



horizons have been the subject of intense scrutiny. The 
de Sitter event horizon considered by Gibbons and Hawk- 
ing as a thermodynamical system |5j is static due to the 
high symmetry of de Sitter space, which admits a timelike 
Killing vector, and plays a role analogous to that of the 
Schwarzschild event horizon among black holes. For more 
general Friedmann-Lemaitre- Robertson- Walker (FLRW) 
spaces, which do not admit such a Killing vector, the 
particle and event horizons are familiar from standard 
cosmology textbooks. However, they do not exist in all 
FLRW spaces and they do not seem suitable for formu- 
lating consistent thermodynamics f [Iol - [l5| and references 
therein). Instead, the FLRW apparent horizon, which al- 
ways exists contrary to the event and particle horizons, 
seems a better candidate. In this paper we reconsider our 
knowledge of this horizon and try to deepen our under- 
standing of it. Specifically, we derive a simple criterion 
for the apparent horizon to be also a trapping horizon and 
we show that the Kodama-Hayward temperature, which 
is based on the Kodama vector playing the role of the 
timelike Killing vector outside the horizon, is positive if 
and only if the apparent horizon is trapping. The causal 
character of this surface is related to this criterion. 

The thermodynamics of the FLRW apparent horizon 
has seen much interest recently, with many authors de- 
riving the temperature of this horizon with the Hamilton- 
Jacobi variant of the Parikh-Wilczek "tunneling" ap- 
proach [l6j]. However, different definitions of surface 
gravity can be applied to this calculation in order to 
define the energy of (scalar) particles, corresponding to 
a background notion of time, and these different pre- 
scriptions provide different notions of temperature. The 
Kodama-Hayward prescription seems to stand out among 
its competitors because the Kodama vector is associ- 
ated with a conserved current even in the absence of 
a timelike Killing vector [rTj . a fact called the "Ko- 
dama miracle" [18| which leads to several interesting re- 
sults. What is more, the Noether charge associated with 
the Kodama vector is the Misner-Sharp-Hernandez mass 
[19L |20( , which is almost universally adopted as the inter- 
nal energy U in horizon thermodynamics. The Misner- 
Sharp-Hernandez mass, defined in spherical symmetry, 
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coincides with the Hawking quasi-local energy 2l| and, 
if we insist in using it in thermodynamics, the use of the 
Kodama-Hayward surface gravity follows naturally (al- 
though this point may be considered debatable by some). 

In the next section, we review background material 
while deriving new formulas useful in the study of appar- 
ent and trapping horizons. Sec. 3 discusses the various 
notions of horizons in FLRW space and their dynamics, 
and elucidates the causal character of the apparent hori- 
zon. The following section raises a question neglected in 
the literature, namely the condition under which the ap- 
parent horizon is also a trapping horizon. The simple cri- 
terion is that the Ricci scalar must be positive. We then 
show that the Kodama-Hayward temperature is positive 
when the apparent horizon is trapping. Sees. 5 and 6 con- 
tain discussions of the thermodynamics of cosmological 
horizons in General Relativity (GR) and in extended the- 
ories of gravity and uses the generalized 2nd law to reject 
as unphysical an exact solution of Brans Dicke theory. 
Sec. 7 contains the conclusions. We follow the notations 
of 22j. The speed of light c, reduced Planck constant 



H, and Boltzmann constant Kb are set equal to unity, 
however they are occasionally restored for better clarity. 

BACKGROUND 

The FLRW line element in comoving coordinates 
(t, r, 9, (p) is 



To transform to the Schwarzschild-like form, one first 
introduces the new time T defined by 



dT 



- (dt + pdR) , 



(3) 



where F is a (generally non-unique) integrating factor 
satisfying 



d_ 

dR 



d_ 

dt 



(4) 



to guarantee that dT is a locally exact differential, while 
/3(i, R) is a function to be determined. Substituting dt — 
FdT — f3dR into the line element, one obtains 



ds 2 



H 2 R 2 



+ 



-11- 



1 — kr 2 
H 2 R 2 



p 2 dT 2 



2 1 



2HRF 



1 — kr 2 



1 — kr 2 

H 2 R 2 
1 — kr 2 

dTdR 



2HR(3 + 1 



1 - kr 2 



dR 2 



Ff3dTdR 



R ' dVl^} 



By choosing 



HR 



1 - H 2 R 2 - kr 2 



(5) 



(6) 



ds 2 = -dt 2 + a 2 it) 



dr 



- 2 dn 2 



1 — kr 2 



(2: 



(i) 



where k is the curvature index, a(t) is the scale factor, 
and dfl 2 2 -, = d9 2 + sin 2 8 dip 2 is the line element on the 
unit 2-sphere. Sometimes different coordinates employ- 
ing the areal radius R(t,r) = a(t)r are useful. Such co- 
ordinate systems include the pseudo-Painleve-Gullstrand 
and Schwarzschild-like coordinates, which we introduce 
here for a general FLRW space. 

Begin from the metric ([T]); using the areal radius R, 
this line element assumes the pseudo-Painleve-Gullstrand 
form 



ds 2 



H 2 R 2 



1 - kR 2 /a 2 
dR 2 



dt 2 



2HR 



1 - kR 2 /a- 



■dtdR 



l-kR 2 /a 2 (2) 



(2) 



where H = a/a is the Hubble parameter and an overdot 
denotes differentiation with respect to the comoving time 
t. We use the word "pseudo" because the coefficient of 
dR 2 is not unity, as required for Painleve-Gullstrand co- 
ordinates [HI and the spacelike surfaces t ^constant are 
not flat (unless k = 0), which is regarded as the essential 
property of Painleve-Gullstrand coordinates 26 1. 



the cross-term proportional to dTdR is eliminated and 
one obtains the FLRW line element in the Schwarzschild- 
like form fH 



ds 2 



H 2 R 2 



l-kR 2 /a 2 
dR 2 



l-kR 2 /a 2 -H 2 R 2 



F 2 dT 2 



(7) 



where F = F(T, R), a, and H are implicit functions of T. 

Horizons in spherical symmetry are discussed in a 
clear and elegant way by Nielsen and Visser in 27 1 
(see also [Hj]). These authors consider the most gen- 
eral spherically symmetric metric with a spherically sym- 
metric spacetime slicing, which assumes the form (in 
Schwarzschild-like coordinates) 



ds 2 



-2<t>(t,R) 



-R 2 d£l 2 2^ 



2M(t, R) 
R 



dt 2 



dR 2 



1 - 



2M(t,H) 



R 



(8) 



where M(t,R) a posteriori turns out to be the Misner- 
Sharp-Hernandez mass 3, 2(|. This form is ultimately 
inspired by the Morris-Thorne wormhole metric [28|, it 
compromises between the latter and the widely used 
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gauge ds 2 = ~A{t,R)dt 2 + B(t, R)dR 2 + R 2 dn 2 2) , and 
is particularly convenient in the study of both static and 
time- varying black holes [27],[2!|. For the metric 0, we 
have 



F(T,R) 
y/l-kR 2 /c 



(9) 



and 



2M kR 2 9 , 8tt 9 

1 - — = 1 - — - H 2 R 2 = 1-— P R 2 
R a 1 3 



i2 „ .. , , , ,10, 

which is consistent with the well known expression 

k \ R 3 4tt 



M = H 



R p 



(11) 



of the Misner-Sharp-Hernandez mass in FLRW space 

Hi. 

In non-spatially flat FLRW spaces, k ^ 0, the quantity 
4tt.R 3 /3 is not the proper volume of a sphere of radius R, 
which is instead 



v proper 



dip d9 / dr V 5 ( 3 ) , (12) 
o Jo Jo 

where g' 3 ' = a \_ s ^i 9 is the determinant of the restric- 
tion of the metric g ab to the 3-surfaces r =constant. 
Therefore, 



4ira 3 (t) 



dr'r' 2 f x 
?LL==±„a 3 (t) d X 'f 2 (x), 
VI - kr' z Jo 

(13) 



where x is the hyperspherical radius and 

sinh x if k < , 



fix) = r = < 



X 



, smx 



if k = , 
if k > , 



(14) 



with x = f = J 



= . Integration gives 



V' 



27ra 3 (t) (rVl + r 2 — sinh 1 r) 



^a 3 (0? 
t 2?ra 3 (t) (si 



if jfe = -1. 
if fc = , 



sin 1 r — r 



) if = +i . 

_ (15) 

However, it turns out that only the "areal volume" 



V 



AnR 3 



(16) 



is used, as a consequence of the use of the Misner-Sharp- 
Hernandez mass, which is identified as the internal energy 
U in the thermodynamics of the apparent horizon. 

The Misner-Sharp-Hernandez mass (fTTj) of a sphere of 
radius R does not depend explicitly on the pressure P 



of the cosmic fluid. Its time derivative, instead, depends 
explicitly on P; consider a sphere of proper radius R = 
R s {t), then, using R = ar and eq. (|3T|) . one has 



M = 4ttR 3 



(17) 



If the sphere is comoving, R s oc a(t), then R s /R s = H 
and 



M = -AttHR 3 P ■ 



(18) 



in this case M depends explicitly on P but not on p. By 
taking the ratio of eqs. (fl"8"|) and (fTTj) one also obtains, in 
GR, 

P 

M + 3H — M = (comoving sphere). (19) 
P 

However, for the thermodynamics of the apparent hori- 
zon (and of the event horizon as well), the horizon is not 
a comoving surface. 

It is now easy to locate the apparent horizon of a gen- 
eral FLRW space. In spherically symmetric spacetimes, 
the apparent horizon (existence, location, dynamics, sur- 
face gravity, etc.) can be studied by using the Misner- 
Sharp-Hernandez mass M [3, [2(| , which coincides with 
the Hawking-Hayward quasi-local mass 0, H3 for these 
spacetimes. The Misner-Sharp-Hernandez mass is only 
defined for spherically symmetric spacetimes. A spheri- 
cally symmetric line element can always be written as 



ds 2 = h ab dx a dx b + R 2 dn 2 {2) , 



(20) 



where a,b = 1, 2 . The Misner-Sharp-Hernandez mass M 
is defined by 



l-^ = V c i?V c i? 
R 



or [87 



M 



R 
~2 



l-h ab V a RV b R) 



(21) 



(22) 



an invariant quantity of the 2-space normal to the 2- 
spheres of symmetry. In a FLRW space, setting g RR = 
(equivalent to /i ab V Q i?Vhi? = in Schwarzschild-like 
coordinates or to Rah = 2M) yields the radius of the 
FLRW apparent horizon 



Ra 



1 



H 



^JH 2 + k/a 2 



(23) 



Since the Misner-Sharp-Hernandez mass is defined quasi- 
locally 21| , this derivation illustrates the quasi-local na- 
ture of the apparent horizon, as opposed to the global 
nature of the event and particle horizons. 
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Eqs. (EJ) and yield 

1- — = 1-^-— (24) 
and the Hamiltonian constraint p5|) then implies that 

(25) 



The Kodama vector is introduced as follows. Using 
the metric decomposition (|20p , let e a b be the volume form 
associated with the 2-metric h a b', then the Kodama vector 
is 



17] 



(26) 



with K 9 = K v = 0. The Kodama vector lies in the 
2-surface orthogonal to the 2-spheres of symmetry and 
K a V a R = e a& V a i?V fe i? = 0. In a static spacetime, 
the Kodama vector is parallel (in general, not equal) 
to the timelike Killing vector. In the region in which 
it is timelike, the Kodama vector defines a class of pre- 
ferred observers with four-velocity u a = K a / ^J\K C K C \. 
It can be proved ([131, see [HI for a simplified proof) 
that the Kodama vector is divergence-free, V a K a = 0, 
which has the consequence that the Kodama energy cur- 
rent J a = G ab K~b is covariantly conserved, V a J a = 0, a 
remarkable property referred to as the "Kodama mira- 
cle" [lij . If the spherically symmetric metric is written 
in the gauge 

ds 2 = -A (t, R) dt 2 + B (t, R) dR 2 + R 2 dn 2 2) , (27) 
then the Kodama vector assumes the simple form (e.g., 

M 



K a 



1 



(28) 



It is shown in [21( that the Noether charge associated 
with the Kodama current is the Misner-Sharp-Hcrnandez 
energy lj|, 2(| of spacetime. The Hay ward proposal for 
the horizon surface gravity in spherical symmetry [l2| is 
based on the Kodama vector. This definition is unique 
because the Kodama vector is unique and KKodama agrees 
with the surface gravity on the horizon of a Reissner- 
Nordstrom black hole, but not with other definitions of 
dynamical surface gravity. The Kodama-Hayward sur- 
face gravity can be written as fl2| 



^Kodama 



1 



2V-h 



WdvR) . (29) 



The components of the Kodama vector in Schwarzschild- 
like coordinates are 



if" = 



yJl-kR 2 /a? 



,0,0,0 



(30) 



and its norm squared is 
K r K c = 



1 _ H 2 R z _ 



R 2 



(31) 

The Kodama vector is timelike (K C K C < 0) if R < Rah, 
null if R = Rah, and spacelike (K C K C > 0) outside the 
apparent horizon R > Rah- 

The components of the Kodama vector in pseudo- 
Painleve-Gullstrand coordinates are 



(32) 



K» = (v/l-/ci? 2 /a 2 ,0,0,o) , 
while in comoving coordinates they are 

K 1 * = (Vl - kr 2 , -HrVl - kr 2 , 0, o) , (33) 



with K C K C = - (1 - kr 2 - a 2 r 2 ) = 1-2M/R (e.g., [33j]). 

In GR, if the FLRW universe is sourced by a perfect 
fluid with energy-momentum tensor 



Tab = (P + p) u a u b + Pg ab , 



(34) 



where p, P, and u a are the energy density, pressure, and 
four-velocity field of the fluid, respectively, one has 



H 



SttG k 



a 
a 



AttG 



(P + 3P) 



(35) 



(36) 



The covariant conservation equation V h T a b = yields 
the energy conservation equation 



p + 3H (P + p) = 



(37) 



which is not independent of eqs. (|35l) and Q36p and can 
be derived from them. Another useful relation following 
from these equations is 



H = -4ttG (P + p) + 4 



(38) 



Let t = denote the Big Bang singularity (in the cases 
in which it is present). All comoving observers whose 
worldlines have u a as tangent are equivalent and, there- 
fore, the following considerations apply to any of them, 
although we refer explicitly to a comoving observer lo- 
cated at r = 0. 



FLRW HORIZONS AND THEIR DYNAMICS 

Two horizons of FLRW space are familiar from stan- 
dard cosmology textbooks: the particle and the event 
horizons The particle horizon [34] at time f is a 



5 



sphere centered on the comoving observer at r 
with radius 



R PH {t) = a{t) 



dt' 
o W)' 



and 



(39) 



The particle horizon contains every particle signal that 
has reached the observer between the time of the Big 
Bang t — and the time t [88| . For particles travel- 
ling radially to the observer at light speed, it is ds = 
and df2(2) = 0. The line element can be written using 
hyperspherical coordinates, 



ds z 



dt 2 +a 2 (t) d X 2 + f 2 (x)dn 2 2) 



(40) 



Along radial null geodesies, d\ = —dt/a and the infinites- 
imal proper radius is a(t)dx)- Integrating between the 
emission of a light signal at \e at time t e and its detec- 
tion at x = at time t, one obtains 



dx 



dt' 



(41) 



and, using Xe — Jq" dx — — f dx, we obtain 



dt' 
W) 



(42) 



The physical (proper) radius R is obtained by multipli- 
cation by the scale factor |89j . 



Re = a(t) 



Now take the limit t e — > + : 



dt' 



(43) 



• if the integral f Q diverges, it is possible for the 
observer at r = to receive all the light signals 
emitted at sufficiently early times from any point 
in the universe. The maximal volume that can be 
causally connected to the observer at time t is infi- 
nite. 

• If the integral J Q -^j^ is finite, the observer at r = 
receives, at time t, only the light signals started 
within the sphere r < J Q j^fry- 

The physical (proper) radius of the particle horizon is 
therefore given by eq. (|3"9l . At a given time t, the particle 
horizon is the boundary between the worldlines that can 
be seen by the observer and those ( "beyond the horizon" ) 
which cannot be seen. This boundary hides events which 
cannot be known by that observer at time t and it evolves 
with time. The particle horizon is the horizon commonly 
studied in inflationary cosmology. 

The particle and event horizons depend on the ob- 
server: contrary to the event horizon of the Schwarzschild 
black hole, different comoving observers in FLRW space 



will see event horizons located at different places. An- 
other difference with respect to a black hole horizon is 
that the observer is located inside the event horizon and 
cannot be reached by signals sent from the outside. 

The cosmological particle horizon is a null surface. 
This statement is obvious from the fact that the event 
horizon is a causal boundary and is generated by the 
null geodesies which barely fail to reach the observer; it 
can also be checked explicitly. Using hyperspherical co- 
ordinates {t,x)i the equation of the particle horizon is 



dt' 

^0 



0. 



The normal to this surface has components 



PH 



= £|ti - 



a 



(44) 



(45) 



and it is straightforward to see that N a N a = 0. 

The particle horizon evolves according to the equation 
(e.g., [13) 



Rph = HRph + 1 - 



(46) 



which is obtained by differentiating eq. (j39J). In an ex- 
panding universe with a particle horizon it is Rph > 0, 
which means that more and more signals emitted between 
the Big Bang and time t reach the observer as time pro- 
gresses. If Rpu(t) does not diverge as t — > t max , then 
there will always be a region unaccessible to the comov- 
ing observers. 

The acceleration of the particle horizon is 

Qj 47T 

Rph = -Rph +H = -—(p + 3P)Rph +H, (47) 
a 3 

Let us turn now our attention to the event horizon. 
Consider all the events which can be seen by the comov- 
ing observer at r — between time t and future infinity 
t = +oo (in a closed universe which recollapses, or in a 
Big Rip universe which ends at a finite time, substitute 
+oo with the time t max corresponding to the maximal 
expansion or the Big Rip, respectively). The comoving 
radius of the region which can be seen by this observer is 



Xeh 



dt' 



(48) 



if this integral diverges as the upper limit of integration 
goes to infinity or to t max , it is said that there is no event 
horizon in this FLRW space and events arbitrarily far 
away can eventually be seen by the observer by waiting 
a sufficiently long time. If the integral converges, there 
is an event horizon: events beyond Teh will never be 
known to the observer 3J] . The physical (proper) radius 
of the event horizon is 



R EH {t) = o(i) 



dt' 
W) 



(49) 
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In short, the event horizon can be said to be the "com- 
plement" of the particle horizon [35|; it is the (proper) 
distance to the most distant event that the observer will 
ever see. Clearly, in order to define the event horizon, 
one must know the entire future history of the universe 
from time t to infinity and the event horizon is defined 
globally, not locally. 

The cosmological event horizon is a null surface. 
Again, the statement follows from the fact that the event 
horizon is a causal boundary. To check explicitly, use the 
equation of the event horizon in comoving coordinates 



d£ 



= 0: 



a(t<) 

the normal to this surface has components 



(50) 



(51) 



and it is easy to see that N a N a = 0. 

The event horizon evolves according to the equation 



Reh = HReh — 1 . 



(52) 



which is obtained by differentiating eq. (|49|). The ac- 
celeration of the event horizon is also straightforward to 
derive, 



Reh — \H + H ) Reh — H . 



(53) 



The event horizon does not exist in every FLRW space. 
To wit, consider a spatially flat FLRW universe sourced 
by a perfect fluid with equation of state P = wp and 
w =const.> —1; if w > —1/3 (i.e., in GR, for a deceler- 
ating universe), there is no event horizon because 



a(t) = a i 3 (™ +1 » 



(54) 



as an estimate of the radius of the event horizon during 
inflation, when the universe is close to a de Sitter space 
4 1| . The Hubble horizon coincides with the apparent 
horizon for spatially flat universes (see eq. (j6"6")l below) 
and with the event horizon of de Sitter space. However, 
this concept does not add to the discussion of the various 
types of FLRW horizons and it seems unnecessary. 

Let us consider now the apparent horizon, which de- 
pends on the spacetime slicing (this feature is illustrated 
by the fact that it is possible to find non-spherical slic- 
ings of the Schwarzschild spacetime without any appar- 
ent horizon 42j,|43j). In a FLRW spacetime, it is natural 
to use a slicing with hypersurfaces of homogeneity and 
isotropy (surfaces of constant comoving time). FLRW 
space is spherically symmetric about every point of space 
and the outgoing and ingoing radial null geodesies have 
tangent fields with comoving components 



„ / VI - kr 2 



(57) 

respectively, as is immediately obtained by setting p c p° = 
for the tangents. There is freedom to rescale a null vec- 
tor by an arbitrary constant (which must be positive if we 
want to keep this vector future-oriented) . The choice (|57| 
implies that l c n c = —2. The more common normaliza- 
tion l c n c = — 1 is obtained by dividing both l a and n a by 
V2. 

The expansions of the null geodesic congruences are 
computed using the equation 



ab 



l a n b 



i a l b 



(-n c l d g cd ) 



VaZfc. 



Computing first 



= 3H 



— y/l — kr 2 , 
ar 



(58) 



(59) 



and the event horizon has radius 



Reh = i 3( ™ +1 » 



3(w + 1) ,_3™+i 
3w + 1 



t 3(^+1) 



(55) 



If w > —1/3, the exponent ^j^m is positive and the in- 
tegral diverges: there is no event horizon in this case. In- 
deed, the existence of cosmological event horizons seems 
to require the violation of the strong energy condition in 
at least some region of spacetime [40] • We can state that 
in GR with a perfect fluid the event horizon exists only 
for accelerated universes with P < —p/3. 

The literature sometimes refers to a "Hubble horizon" 
of FLRW space with radius 



Rh = 



1 

H 



(56) 



This quantity only provides the order of magnitude of 
the radius of curvature of a FLRW space and is used 



V c n c = ZH v 7 ! ~ kr 2 , 

ar 




a?r 2 sin 2 i 



9cdl c n d = -2 , 



VI - kr 2 ' 

r^ = o, r§! = r^ = m cl , r c n 

the result is [90j 



ar 



R 



(60) 



krS ci +aa5 
1 — kr 2 



cO 



% = 2(ar + VT=W) =2 f g+ l^*^ 



(61) 



2 (ar - VI - kr 2 ) / 1 / kR 2 

ar = 2 r ~:r v 1 - — 



(62) 
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Following [44j, the apparent horizon is a surface defined 
by the conditions on the time slicings 



>0, 

= 0, 



and is located at 



rAH 



or 



Rah {t) 



Va 2 + k 



y/H 2 + k/a 2 



(63) 
(64) 



(65) 



(66) 



in terms of the proper radius R = ar. The apparent 
horizon is defined locally using null geodesic congruences 
and their expansions, and there is no reference to the 
global causal structure. 

Looking at eqs. (I5T1) and ([B"2")) . or at their product 



R 2 \R\h 



(67) 



it is clear that when R > Rah it is 9i > and 9 n > 0, 
while the region < R < Rah has 9i > and 9 n < 
(radial null rays coming from the region outside the 
horizon will not cross it and reach the observer). 

For a spatially flat universe, the radius of the appar- 
ent horizon Rah coincides with the Hubble radius H , 
while for a positively curved (fc > 0) universe Rah is 
smaller than the Hubble radius, and it is larger for an 
open (k < 0) universe. In GR, the Hamiltonian con- 
straint (|35j) guarantees that the argument of the square 
root in eq. (|66|) is positive for positive densities p. The 
apparent horizon exists in all FLRW spaces. 

In general, the apparent horizon is not a null surface, 
contrary to the event and particle horizons. The equation 
of the apparent horizon in comoving coordinates is 



F{t,r) = a{t)r 



y/H 2 + k/c 



0. 



(68) 



The normal has components 



tyi = V^\ah 



h(h -k/a< 
{H 2 + k/a 2 ) 3/2 



Sao + adui 



AH 



= HRah 



H ~~^2l R ' AB 



5uq + aS l 



Mi 



HRah ~ <^o + a "W • 
a 



(69) 



In GR with a perfect fluid with equation of state P — wp, 
eqs. dH]), (HI, and flU yield 



N„ 



(3w + 1) 



HRahS, 



[i0 



(70) 



The norm squared of the normal is 



N a N a = \-kr\ H - 



H 2 H + H' 



{H 2 + k/a 2 f 



H 2 R 2 



AH 



R 



4 

AH 



3H 2 R\ H {p + p){p _ 



4p 2 



H 2 R\h (1 - ^H^Rah) ) 



3P) 



2 rr4 D 4 



(71) 



(72) 



where q = ~aa/a 2 is the deceleration parameter. The 
horizon is null if and only if P = — p or P = p/3. In 
GR with a perfect fluid, the Hamiltonian constraint ([33)) 
yields 



2 D 2 



H l R 



AH 



k 



i 2 H 2 



{8ttG 



Pc 

p 



(73) 

where p c = is the critical density and = p/p c is 
the density parameter, and one obtains 



N a N a = ~(w + l)(3w- 1) H 2 R 2 AH = 
Eq. ([74]) establishes that: 



(74) 



• if —1 < w < 1/3, then N C N C > and the apparent 
horizon is timelike. For a k — universe in Ein- 
stein's theory this condition corresponds to H < 0. 

• If w = -1 or w = 1/3, then N C N C = and the 
apparent horizon is null (de Sitter space, which has 
H = and q = — 1 , falls into this category but it is 
not the only space with these properties). 

• If w < — 1 or w > 1/3, then iVW c < 0, the nor- 
mal is timelike, and the apparent horizon is space- 
like. In Einstein's theory with k — and a perfect 
fluid as the source, w < — 1 corresponds to H > 
("superacceleration"). This is the case of Big Rip 
universes and of a phantom fluid which violates the 
weak energy condition. 

The black hole dynamical horizons considered in the 
literature are usually required to be spacelike [9]. How- 
ever, cosmological horizons can be timelike. In GR, 
the radius of the apparent horizon can be written as 



in terms of the density parameter 



Rah = (Vn \H\ 
VI by using eq. ([73]) . 

The a ppa rent horizon evolves according to the equa- 



tion 



31,11-1481 



k 



Rah = HRah \^-H)= ^HR ah (P + p) , (75) 
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as is easy to check by differentiating eq. (I66t with respect 
to t. fn GR with a perfect fluid as a source, the only way 
to obtain a stationary apparent horizon is when P = 
—p. For de Sitter space, eq. (1751) reduces to Rah = 0, 
consistent with Rh = H^ 1 and H =const. For non- 
spatially flat universes, the equation of state P = —p 
produces other solutions. For example, for k = — 1 and a 
cosmological constant A > as the only source of gravity, 
the scale factor 



apparent horizons 



(76) 



is a solution of the Einstein-Friedmann equations. The 
radius of the event horizon has the time dependence 



R E H{t) =Jj sinh ( sj ~t 



In 




(77) 

The apparent horizon, instead, has constant radius 
Rah = 

As another example consider, for k = +1 and cosmo- 
logical constant A > 0, the scale factor 



a(t) 

the event horizon has radius 
REn{t) = 




(78) 



i C0Sh (Vt* 



7T 



tan 




(79) 



where n = 0,±1,±2, ... The multiple possible values of 
n correspond to the infinite possible branches which one 
can consider when inverting the tangent function, and to 
the fact that in a closed universe light rays can travel 
multiple times around the universe. In this situation it is 
pro blematic to regard the event horizon as a true horizon 
[491 ] . The apparent horizon has constant radius Rah — 

r, according to the fact that pa + -Pa = in eq. (|75|) 

In a k = FLRW universe with a perfect fluid and con- 
stant equation of state P = wp and —1 < w < —1/3 (ac- 
celerating but not superaccelerating universe), the event 
horizon is always outside the apparent horizon and is, 
therefore, unobservable 14, 50| . 

Let us summarize the dynamical evolution of the 
FLRW horizons and compare their evolutionary laws. 
The first question to ask is whether these horizons are 
comoving: they almost never are. The difference be- 
tween the expansion rate of a horizon R/ R and that of 
the expanding matter H is, for the particle, event, and 



Rph 
Rph 

Reh 
Reh 

Rah 
Rah 



H 



1 



-H = - 



Rph 
1 



Reh 



-H = H 



H 2 + -! 

a 

(Sw + l)H 



- 1 



H R 



(80) 
(81) 

2 

AH 

(82) 



respectively. Taking into consideration only expanding 
FLRW universes (Pt > 0), when it exists the particle 
horizon always expands faster than comoving. The event 
horizon (which only exists for accelerated universes) al- 
ways expands slower than comoving. The apparent hori- 
zon expands faster than comoving for decelerated uni- 
verses (a < 0); slower than comoving for accelerated 
universes (a > 0); and comoving for coasting universes 
(a(t) oc i). 

An even simpler way of looking at the evolution is by 
using the comoving radius of the horizon: if this radius 
is constant, then the horizon is comoving. We have, 



r PH = - > , 
a 



Teh = < , 

a 



r AH 



aa 



(a 2 + k) 



3/2 



(83) 
(84) 
(85) 



respectively. The causal character and the dynamics of 
the various FLRW horizons are summarized in Tables U 
andHH 



TRAPPING HORIZON OF FLRW SPACE 

Let us now ask the question: When is the FLRW ap- 
parent horizon also a trapping horizon? According to 
Hayward's definition, when Ci9 n > 0, which gives the 
coordinate- (but not slicing)-invariant criterion 



H a 



>0, 



(86) 



where R a a is the Ricci scalar of FLRW space, and Ci 
is the Lie derivative along l a . In fact, using eqs. J3J 



9 



and (|62|) . we have 

£i6n — ^ a V a #„ = l a d a t 



2 dt+ y^^ dr ) h-^E^ 



— I HE 2 + HR\ I 



fci? 2 



At the apparent horizon R = Rah it is 

ClOn \ A H 



25+ 2ff 2 



H 



H 



+ 



2 {H 2 +k/a 11 



H 2 + £ ^jm + k/a 2 

k \ i? Q a 



This result is independent of the field equations. If we 
assume Einstein's theory and a perfect fluid as the sole 
source of gravity, we obtain 



4*nU= ^(P"3P) 



(87) 



and, therefore, 



i/ie apparent horizon is also a trapping horizon iff R a a > 
(equivalent to P < p/3 in GR with a perfect fluid) . 

Note that in a radiation-dominated universe, which is 
decelerated, the event horizon does not exist. 



THERMODYNAMICS OF COSMOLOGICAL 
HORIZONS IN GR 

Originally developed for static or stationary event hori- 
zons, black hole thermodynamics has now been extended 
to apparent, trapping, isolated, dynamical, and slowly 
evolving horizons [7H9J . Similarly, the cosmological ther- 
modynamics associated with cosmological horizons has 
been extended from the static de Sitter event horizon [B| 
to FLRW dynamical apparent horizons. 

The thermodynamic formulae valid for the de Sitter 
event (and apparent) horizon are generalized to the non- 
static apparent horizon of FLRW space. The apparent 
horizon is argued to be a causal horizon associated with 
gravitational temperature, entropy and surface gravity 
in dynamical spacetimes ( T3~ 15 1 and references therein) 
and these arguments apply also to cosmological horizons. 
That thermodynamics are ill-defined for the event hori- 



zon of FLRW space was argued in [14 . |49H51|- The 
Hawking radiation of the FLRW apparent horizon was 



computed in 52|, [53 1. The authors of 1231 125| rederived 
it using the Hamilton- Jacobi method [27l l54l IHBj in the 
Parikh-Wilczek approach originally developed for black 
hole horizons (16|. In this context, the particle emission 
rate in the WKB approximation is the tunneling proba- 
bility for the classically forbidden trajectories from inside 
to outside the horizon, 



exp 



2Im(I)\ 



exp 



Hid \ 

KbTJ 



where / is the Euclideanized action with imaginary part 
Im(/) , u) is the angular frequency of the radiated quanta 
(taken, for simplicity, to be those of a massless scalar 
field, which is the simplest field to perform Hawking effect 
calculations), and the Hawking temperature is read off 
the expression of the Boltzmann factor, KbT = - figr^ ■ 
The particle energy hw is defined in an invariant way as 
u = —K a V a I, where K a is the Kodama vector, and the 
action / satisfies the Hamilton- Jacobi equation 



h ab \7 a IV b I = 0. 



(89) 



Although the definition of energy is coordinate-invariant, 
it depends on the choice of time, here defined as the Ko- 
dama time. 

A review of the thermodynamical properties of the 
FLRW apparent horizon, as well as the computation of 
the Kodama vector, Kodama-Hayward surface gravity, 
and Hawking tem per ature in various coordinate systems 
are given in Ref. [331 ] . The Kodama-Hayward tempera- 
ture of the FLRW apparent horizon is given by 



K B T 



Rah[H z + ^ 



k 
■2tt- 



2tt 



tt; — ) RAnR a a — ( — 

247TC / V c 



Ra 



H 



(p-3P) . 
(90) 



The expression of the temperature depends on the choice 
of surface gravity k since T = \k\/2tt in geometrized units 
and there are several inequivalent prescriptions for this 
quantity (see 5|1 57 1 for reviews). The choice of k giving 
the temperature reported here is the Kodama-Hayward 
prescription ([29]) [33j. In fact, this equation yields 



^Kodama — y— ^2iJ 2 + H + —^J — R" a , 

(91) 

as can be quickly assessed by using comoving coordi- 
nates and the decomposition (|20|) of the metric, where 

h ab =diagf— 1, 1 ®l r i \ The entropy of the FLRW ap- 
parent horizon is 



Sah = 













V SG } 



, (92) 
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where 



A AH = 4ttR 2 ah = 



-in 



H 2 + k/a? 



(93) 



is the area of the event horizon. The Hamiltonian con- 
straint (f3"5l gives 



Sah = 



S A 



9H 



H 



(P + p) 



(94) 



In an expanding universe the apparent horizon entropy 
increases if P + p > 0, stays constant if P = —p, and 
decreases if the weak energy condition is violated, P < 
-P- 

It seems to have gone unnoticed in the literature that 
the horizon temperature is positive if and only if the Ricci 
scalar is, which is equivalent to equations of state satis- 
fying P < p/3 for a perfect fluid in Einstein's theory. 
This is the condition for the apparent horizon to be also 
a trapping horizon. A "cold horizon" with T — is ob- 
tained for vanishing Ricci scalar but the entropy is pos- 
itive for such an horizon, a situation analogous to that 
of extremal black hole horizons in GR. Note also that, 
if the weak energy condition (which implies P + p > 0) 
is assumed, the boundary P = p/3 between positive and 
negative Kodama-Hayward temperatures corresponds to 
the boundary between timelike and spacelike character 
of the apparent horizon. It is not obvious a priori that a 
null apparent horizon, obtained for R a a = (P = p/3 for 
a perfect fluid in GR) , should occur when the universe is 
filled with conformal matter. 

The natural choice of surface gravity seems to be 
that of Kodama-Hayward, which produces the appar- 
ent horizon temperature (|90[) . The apparent horizon en- 
tropy is Sah — : 4 C = ' k P\h (^ ms can be obtained 
using Wald's Noether charge method — see the discus- 
sion of the next section), and the internal energy U 
should be identified with the Misner- Sharp-Hernandez 

mass Mah — 47tF ^ ah p contained inside the apparent 
horizon. The factor 4irR AH /3 is not the proper volume 
of a sphere of proper (areal) radius Rah unless the uni- 
verse has flat spatial sections. It is the use of the Misner- 
Sharp-Hernandez mass which points us to use the areal 
volume Vah = AlrR ^ AH instead of the proper volume when 
discussing thermodynamics (failing to do so would jeop- 
ardize the possibility of writing the 1st law consistently). 
However, even with this caveat, the 1st law does not as- 
sume the form 



TahSah = Mah + PVah 



(95) 



that one might expect. Let us review now the laws of 
thermodynamics for cosmological horizons. 

0th law. The temperature (or, equivalently, the surface 
gravity) is constant on the horizon. This law ensures that 
all points of the horizon are at the same temperature, or 



that there is no temperature gradient on it. The 0th law 
is a rather trivial consequence of spherical symmetry. 

1st law. The 1st law of thermodynamics for apparent 
horizons is more complicated than (|95[) and was given in 
Refs. Ill, [l2[ under the name of "unified 1st law". While 
using the Misner-Sharp- Hernandez mass Mah as internal 
energy, the Kodama-Hayward horizon temperature (|90|) , 
and the areal volume, one introduces further quantities 
as follows [III Ell ■ Decompose the metric as in eq. (|20p ; 
then the work density is 



w = -^T ab h ab 



ip a = T a b V b R + wV a R 



(96) 



(97) 



is the energy flux across the apparent horizon, when 
computed on this hypersurface. The quantity Aah^Po, 
is called the energy supply vector. The quantity 

j a = ip a + wK a (98) 

is a divergence-free energy-momentum vector which can 
be used in lieu of ip a . The Einstein equations then give 

HQ 



M = kR 2 + 4itR 3 w . 



V a M = Aj a . 



(99) 



(100) 



The last equation is rewritten as [llj, [12 1 

Mp a = V a M - wV a V A H (101) 

("unified 1st law"). The energy supply vector is then 
written as 



(102) 



Along the apparent horizon, it is Mah = Rah/% and 



AaH^Po, = 7T- VcSah = TahVuSah 
Zir 



(103) 



This equation is interpreted by saying that the energy 
supply across the apparent horizon AAH^a is the "heat" 
Tah^ ' a^AH gained. Writing the energy supply explicitly 
gives 



Tah^clSah = VqMa^ - vNoVah 



(104) 



and —wVoYah is a work term. The "heat" entering the 
apparent horizon goes into changing the internal energy 
Mah and performing work due to the change in size of 
this horizon. 

Let us compute now the time component of eq. (| 104|) 
in comoving coordinates for a FLRW space sourced by a 
perfect fluid in GR. We have 



W 



[(P + p) u a u b + Pg ab ] h ab = P -^- , (105) 
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Vah = ZHVah \l--R AH \ = 



9HV AH 
2p 



(P + P) 



(106) 



Mah = j t (Vahp) = i -^f L (P + p), (107) 



and 



S AH = 2ttR ah Rah = H{P + p) (108) 



so that 

TahSah 



HVah f l-*R AH )(p-3P) 



3HV AH 
4p 



(P + p) (p-3P) 



Therefore, it is [ll|, [12 1 



TahSah — Mah + - — ^r— Vah ■ 



(109) 



(110) 



In the infinitesimal interval of comoving time dt the 
changes in the thermodynamical quantities are related 
by 



T A HdS A H = dM AH +dW AH , dW AH 



(P-P) , T/ 
2 dVAH ■ 

(111) 

The coefficient of (IV ah > i.e., — w = (P — p) /2 equals 
the pressure P (the naively expected coefficient) only 
if P = —p (which includes de Sitter space in which 
dM A H , dV A H , and dS A H all vanish) . The fact that the 
coefficient appearing in the work term is not simply P 
can be understood as a consequence of the fact that the 
apparent horizon is not comoving. For a comoving sphere 
of radius R s it is R s /R s = H and V s — 3HV S , while 

M s = V sP + V s p = 3HV a p - 3HV S (P + p) = -3HV S P , 

(112) 

hence M s + PV S = 0. Indeed, the covariant conservation 
equation (|37[) is often presented as the 1st law of ther- 
modynamics for a comoving volume V. Because of spa- 
tial homogeneity and isotropy there can be no preferred 
directions and physical spatial vectors in FLRW space, 
therefore the heat flux through a comoving volume must 
be zero. In fact, consider a comoving volume V c (which, 
by definition, is constant in time) and the correspond- 
ing proper volume at time t, V = a?(t)V c . Multiplying 
eq. ([571) by V one obtains 



Vp + V (P + P ) 



PV = 0. 



(113) 



By interpreting U = pV as the total internal energy of 
matter in V , one obtains the relation between variations 
in the time dt 



and the 1st law (with work term coefficient P) then gives 
TdS = 0, which is consistent with the above-mentioned 
absence of entropy flux vectors and with the well known 
fact that, in curved space, there is no entropy generation 
in a perfect fluid (the entropy along fluid lines remains 
constant and there is no exchange of entropy between 
neighbouring fluid lines [HI]). Indeed, eq. (fT9|) for the 
evolution of the Misner-Sharp-Hernandez mass contained 
in a comoving sphere reduces to M + PV = or p + 
3H (P + p) = 0. However, for a non-comoving volume, 
the work term is more complicated than PdV . 

Attempts to write the 1st law for the event-instead of 
the apparent, horizon lead to inconsistencies [14l . 
This fact supports the belief that it is the apparent hori- 
zon which is the relevant quantity in the thermodynamics 
of cosmological horizons. 

(Generalized) 2nd law. A second law of thermodynam- 
ics for the event horizon of de Sitter space wasgiven 
already in the original Gibbons-Hawking paper Q and 
re-proposed in [59| . Davies [4!| has considered the event 
horizon of FLRW space and, for GR with a perfect fluid 
as the source, has proved the following theorem: if the 
cosmological fluid satisfies P + p > and a(t) — > +oo 
as t — > +oo, then the area of the event horizon is non- 
decreasing. The entropy of the event horizon is taken 

to be Seh = { ^na ) ^P"' where A E h is its area. The 
validity of the generalized 2nd law for certain radiation- 
filled universes was established in 

Due to the difficulties with the event horizon one is led 
to consider the apparent horizon instead. Then, eq. (|108[) 
tells us that, in an expanding universe in Einstein's the- 
ory with perfect fluid, the apparent horizon area increases 
except for the quantum vacuum equation of state P = —p 
(for which Sah stays constant) and for phantom fluids 
with P < —p, in which case Sah decreases, adding an- 
other element of weirdness to the behaviour of phantom 
matter (e.g., [(32| and references therein). 

The generalized 2nd law states that the total entropy of 
matter and of the horizon Stotai = S ma tter + Sah cannot 
decrease in any physical process, 



8S = SS n 



SSah > . 



(115) 



dU + PdV = , 



(114) 



(We refer here to the apparent horizon, but several au- 
thors refer instead to the event or particle horizons. The 
apparent horizon is more appropriate since it is a quasi- 
locally defined quantity.) 



ENTROPY OF THE APPARENT HORIZON AND 
SCALAR-TENSOR GRAVITY 

Black hole thermodynamics has been studied in scalar- 
tensor and other theories of gravity (see [63] for a sum- 
mary and a list of references). Numerical studies show 
that in the intermediate stages of collapse of dust to a 
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black hole in Brans-Dicke gravity, the horizon area de- 
creases and the apparent horizon is located outside the 
event horizon [64j . Horizon entropy in Brans-Dicke grav- 
ity was analyzed by Kang [65], who pointed out that 
black hole entropy in this theory is not simply one quar- 
ter of the horizon area, but rather 



~BH 



(116) 



where <j) is the Brans-Dicke scalar field (assumed to be 
constant on the horizon) and 

i s the determinant of 

(2) 

the restriction g^J = g^ | E of the metric g^ to the 
horizon X. Naively, this expression can be understood 
by replacing the Newton constant G with the effective 
gravitational coupling 



G 



eff 



(117) 



of Brans-Dicke theory; then, the quantity Sbh is non- 
decreasing during black hole collapse [65| . Eq. (11161) has 
now been derived using various procedures 

As done in (65| , consider the Einstein frame representa- 
tion of Brans-Dicke theory given by the conformal rescal- 
ing of the metric 



9^ 



9^ 



G<f>, 



accompanied by the scalar field redefinition 
(f> given by 



2uj 



16ttG 



(118) 
!> with 

(119) 



The Brans-Dicke action [6J 



Ibd = I d A x 



16tt 



<j>R-- gTV^Vut - V{4>) + £ (m) 

(120) 

(where £( m ) is the matter Lagrangian density) is mapped 
to its Einstein frame form 



Ibd 



d x y/ —g 

£(m) 



I? 



16ttG 



(Gcf>y 



(121) 



where a tilde denotes Einstein frame quantitites and 



U 



(gH4>)) 2 



(122) 



where = </>(</>)• I n the Einstein frame the gravitational 
coupling is constant but matter couples explicitly to the 
scalar field and massive test particles following geodesies 
of the Jordan frame g^ do no longer follow geodesies of 



g^ v in the Einstein frame. Null geodesies are not changed 
by the conformal transformation. A cosmological or 
black hole event horizon, being a null surface, is also 
unchanged. The area of an event horizon is not, and the 
change in the entropy formula Sbh = jfj — > t<t~jj = 
is merely the change of the horizon area due to the con- 
formal rescaling of g^. In fact, gffi — Q 2 gffl and, since 
the event horizon is not changed, the Einstein frame area 
is 



A = 



Gcj)A (123) 



assuming that the scalar field is constant on the horizon 
(if this is not true the zeroth law of black hole thermody- 
namics will not be satisfied). Therefore, the entropy-area 
relation for the event horizon Sbh = A /AG still holds 
in the Einstein frame. This is expected on dimensional 
grounds since, in vacuo, the theory reduces to GR with 
varying units of length l u ~ fll u , time t u ~ fH u , and 



mass m„ = fl 



(where t u ,l u , and m u are the con- 



stant units of time, length, and mass in the Jordan frame, 
respectively). Derived units vary accordingly [70j. An 
area must scale as A ~ SI 2 = G<f> and, in units in which 
c = h = 1 the entropy is dimensionless and, therefore, is 
not rescaled. As a result, the Jordan frame and Einstein 
frame entropies coincide 



65] 



The equality between black hole entropies in the Jor- 
dan and Einstein frames is indeed extended to all theories 



with action J d 4 



which admit 



an Einstein frame representation [71[. 

As a byproduct of this observation, the Jordan and the 
Einstein frames are physically equivalent with respect to 
the entropy of the event horizon. A debate on whether 
these two frames are physically equivalent seems to flare 
up now and again; it is pretty well established that, at 
the classical level, the two frames are simply different 



72 



Poten- 



representations of the same physics 
tial problems arise from the representation-dependence 
of fundamental properties of theories of gravity (includ- 
ing the Equivalence Principle), but this is not an ar- 
gumente against the equivalence of the two conformal 
frames: rather, it means that such fundamental proper- 
ties should, ideally, be reformulated in a representation- 
independent way ( 7J] and references therein). We will 
not address this problem here. 

The classical equivalence is expected to break down 
at the quantum level; in fact, already in the absence of 
gravity, the quantization of canonically related Hamil- 
tonians produces inequivalent energy spectra and eigen- 
functions (75|. However, it is not clear that this happens 
at the semiclassical level for conformally related frames 
[72] ]. Black hole thermodynamics is not purely classical 
(the Planck constant h appears in the expressions of the 
entropy and temperature of black hole and cosmological 
horizons). It is not insignificant that the physical equiv- 
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alence between conformal frames holds for the (semiclas- 
sical) entropy of event horizons. 

Contrary to event horizons, the location of apparent 
horizons (which, in general, are not null surfaces), is 
changed by conformal transformations. The problem of 
relating Einstein frame apparent horizons to their Jor- 
dan frame counterparts, raised in 76j, has been solved in 
& 

At this point, one may object that there was a logical 
gap in our previous discussion: following common prac- 
tice, we took the formula S = A /AG for the (black hole 
or cosmological) event horizon and we used it for the ap- 
parent horizon. Let us consider stationary black hole or 
cosmological event horizons first. The area formula can 
be derived (in GR or in other theories of gra vity) by using 
Wald's Noether charge method 6^-68, ZS, 23 or other 
methods 0, [f36j]. As a result, the usual entropy-area re- 
lation remains valid provided that the gravitational cou- 
pling G is replaced by the corresponding effective gravi- 
tational coupling G e f f of the theory, the identification of 
which follows from the inspection of the action or of the 
field equations rewritten in the form of effective Einstein 
equations. More rigorously, in [8(| G e // is identified by 
using the matrix of coefficients of the kinetic terms for 
metric perturbations [8oT ] . The metric perturbations con- 
tributing to the Noether charge in Wald's formula are 
identified with specific metric perturbation polarizations 
associated with fluctuations of the area density on the 
bifurcation surface E of the horizon (in D spacetime di- 
mensions, this is the (D — 2)-dimensional spacelike cross- 
section of a Killing horizon on which the Killing field 
vanishes, and coincides with the intersection of the two 
null hypersurfaces comprising this horizon). The horizon 
entropy is 



lated to be feC 



.4 



~>BH 



for a theory described by the action 



(124) 



1= J d A X y/^g C{g llv ,R a i3pX,V ' <rR a p P xA^ a4>, ••■) , 

(125) 

where is a gravitational scalar field. The Noether 
charge is 



S 



2tt / d 2 x 




(126) 



where e pCT is the (antisymmetric) binormal vector to the 
bifurcation surface E (which satisfies V^Xj, = e Mt , on the 
bifurcation surface E, where x M is the Killing field van- 
ishing on the horizon) and is normalized to e ab i a b = — 2. 
The subscript (0) denotes the fact that the quantity in 
brackets is evaluated on solutions of the equations of mo- 
tion. The effective gravitational coupling is then calcu- 



G 



eff 



= -2tt 



SR 



pvper / (o) 



(127) 



For dynamical black holes, there is no timelike Killing 
vector to provide a bifurcate Killing horizon. However, 
it was shown in [111 ] that, for apparent horizons, the Ko- 
dama vector can replace the Killing vector in Wald's en- 
tropy formula, and the result is one quarter of the area 
in GR (or the corresponding generalization in theories 
of the form (|125[1 ). We do not repeat the calculation 
here, but we simply note that the same calculation ap- 
plies to cosmological apparent horizons as well, a point 
that seems to not have been noted in the literature on 
cosmological horizons. 

Finally, let us see how thermodynamics can restrict 
the range of physical solutions of a theory of gravity. 
In Brans-Dicke cosmology, consider the exact solution 
representing a spatially flat universe with parameter u) = 
—4/3, and no matter [81| 



a(t) = ao cxp (Hi) , 



<f>(t) = o exp(-3fft) , 



(128) 



(129) 



where ao, 4>a, and H are positive constants. In GR, de 
Sitter spaces are obtained with constant scalar fields but 
this is not always the case in scalar-tensor gravity. Since 
the entropy of the apparent/event horizon in this case 
is S = (J>Ah/A = <f)oH~ 2 exp(—3Ht), it is always de- 
creasing. The scalar field plays the role of an effective 
fluid with density and pressure given by the equations of 
Brans-Dicke cosmology which, in the spatially flat case 
and for vacuum and a free Brans-Dicke scalar, reduce to 




3H<P = . 



(130) 



(131) 



(132) 



In our case, the energy density and pressure of the effec- 



tive fluid are pu) — - 
this effective fluid is 



5 (« 



H . The entropy density of 



P(<t>) + P{4>) 
T 



= 0; 



(133) 



therefore, the horizon associated with this solution vio- 
lates the generalized 2nd law and should be regarded as 
unphysical. More generally, the use of entropic considera- 
tions to select extended theories of gravity, was suggested 
in 0. 
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CONCLUSIONS 

The apparent horizon suffers from the dependence on 
the spacetime slicing. In FLRW space, it would be unnat- 
ural to choose a slicing unrelated to the hypersurfaces of 
spatial homogeneity and isotropy and constant comov- 
ing time, because the latter identify physical comoving 
observers who see the cosmic microwave background ho- 
mogeneous and isotropic around them (apart from small 
temperature anisotropies of the order 5 TO -5 ). The prob- 
lem of the slicing dependence, therefore, does not seem so 
pressing in FLRW spaces, however it is not completely 
eliminated. Nevertheless, apparent horizons seem bet- 
ter candidates for thermodynamical considerations than 
event or particle horizons. 

Similar to dynamical black hole horizons, the thermo- 
dynamics of FLRW cosmological horizons is not com- 
pletely free of problems: the choice of surface gravity 
determines the horizon temperature and there are many 
inequivalent proposals for surface gravity. A natural 
choice of internal energy contained within the apparent 
(or even event) horizon is given by the Misner-Sharp- 
Hernandez mass, and then it is natural to choose the 
Kodama time and Kodama-Hayward surface gravity be- 
cause the Misner-Sharp-Hernandez mass is intimately as- 
sociated with the Kodama vector as a Noether charge. 
However, doing so, produces a Kodama-Hayward tem- 
perature which is negative for GR universes with stiff 
equations of state P > p/3, and it is not obvious that 
one should give up these universes as unphysical. Differ- 
ent choices of temperature would produce different forms 
of the 1st law, corresponding to different coefficients for 
the work term dW appearing there. What is certain, 
though, is that the apparent (and also the event and 
particle) horizons are not comoving and one should not 
necessarily expect a simple PdV term to appear. Perhaps 
other choices of quasi-local energy can produce consistent 
forms of the 1st law and be applicable to a larger variety 
of universes; one could turn the argument around and use 
cosmology to help selecting the "correct" surface gravity 
also for dynamical black hole horizons (this possibility 
will be the subject of a separate publication). 

We have elucidated the causal character of the appar- 
ent horizon and given a simple criterion for the apparent 
FLRW horizon to be trapping. This criterion coincides 
with the one for the Kodama-Hayward temperature to 
be positive-definite, and the threshold between trapping 
and untrapping horizons is a FLRW universe filled with 
conformal matter which, if the weak energy condition is 
assumed, also marks the transition between timelike and 
spacelike nature of the apparent horizon. At the moment, 
we are unable to offer a simple and consistent physical 
interpretation of this fact and we will refrain from doing 
so. We have also considered the extension of the thermo- 
dynamics of cosmological horizons to alternative theories 



of gravity and, within Brans-Dicke theory, we have seen 
how thermodynamical considerations can help judging 
how physical a certain solution can be. 

Overall, it appears that the thermodynamics of cos- 
mological apparent horizons exhibits features which are 
not yet fully understood. Due to the extremely simpli- 
fied nature of FLRW spacetime, understanding these as- 
pects for cosmological apparent horizons should be more 
fruitful and rapid than understanding the corresponding 
aspects of black hole dynamical horizons. 

This work is supported by the Natural Sciences and 
Engineering Research Council of Canada (NSERC). 
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Horizon 


causal character 


Event horizon 


null 


Particle horizon 


null 


Apparent horizon 


timelike if — p < P < p/3, 
null if P = —p or p/3, 
spacelike if P < — p or P > p/3 


de Sitter horizon 


null 



TABLE I: Causal character of the FLRW cosmological horizons. 
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Horizon 


location 


velocity 


acceleration 


Event horizon 


Reh = a(t) J t + °° 


Reh = HReh — 1 


Reh = ^Reh — H 


Particle horizon 




Rph = HRph + 1 


Rph = ^Rph + H 


Apparent horizon 


Rah yj^W 


D _ 17 D3 / fc _ i \ 

rtAH — nix AH y-^j ±j 
= 4nHR AH (P + p) 


Rah = Rah 
r . / ■ ^\ 

-H + 3H 2 R 2 AH 


de Sitter horizon 


-RdS = H^ 1 


R dS = 


ifcs = 



TABLE II: FLRW cosmological horizons and their dynamical behaviour. 



